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Abstract 

Local magnetic properties around a nonmagnetic impurity are investigated 
in high-Tc cuprate superconductors. We consider a model two-dimensional 
dx2_y2-wme superconductor with strong antiferromagnetic (AF) spin fluctua- 
tions and calculate the spin susceptibility around the impurity. We show that 
the uniform susceptibility (x(0)), which is usually suppressed in spin-singlet 
superconductivity, is enhanced at low temperatures around the impurity when 
the impurity scattering is in the unitarity limit. This anomalous behavior is 
more remarkable when the impurity potential is extended spatially and the 
AF spin fluctuations are stronger. As the origin of this phenomenon, we point 
out the importance of (1) a low-energy bound state around the impurity and 
(2) a mode-mode coupling effect which connects the AF spin fluctuations with 
x(0). Our results can explain the anomalous temperature dependence of the 
Knight shift on Li which was recently observed in YBa2(Cui_a;Lia;)306+3/. 
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I. INTRODUCTION 



Recently, Bobroff et al. observed an anomalous behaviour of the uniform susceptibility 
(x(0)) in YBa2(Cui_2:Li^.)306+j,: ^ Usually x(0) is suppressed in singlet superconductivity 
because of the singlet Cooper pairs. However, their Knight shift measurement clarified that 
x(0) on Li increases with decreasing temperature in the superconducting state. In addition, 
this phenomenon was shown to be more remarkable in the under-doped regime than in the 
over-doped one. Since antiferromagnetic (AF) spin fluctuations are considered stronger in 
the former regime, this result implies the importance of the AF spin fluctuations in this 
phenomenon. 

In our previous papers, we have clarified two nonmagnetic impurity effects on the 
magnetic properties of high-Tc cuprates: 

(A) In d-wave superconductors, when the impurity scattering is in the unitarity limit, a low- 
energy bound state is formed around the impurity. Then, at low temperatures, the 
low-energy AF spin fluctuations are enhanced by the density of states (DOS) originat- 
ing from this bound state. [0] This impurity effect can explain the enhancement of the 
low-energy AF spin fluctuations in Zn-doped YBa2Cu306+x- It also agrees with 
the increase of the nuclear spin-lattice relaxation rate on Cu in YBa2(Cui_2:Zn^)408 
far below the superconducting transition temperature (Tc). [P|-[T^ 



(B) In the normal state, when the local density of states around the Fermi level is enhanced 
near an impurity, the AF spin fluctuations are also enhanced locally. This enhance- 
ment affects x(0) through a mode-coupling effect. Then x(0) shows a Curie- Weiss-like 
temperature dependence around the impurity as observed in YBa2(Cui_sM2:)306+y 
(M=Zn, Li) in the normal state. [14,15|| Recently, this mechanism was discussed also 



by Bulut. [1T6|,|T3 

Taking into account these effects, we can draw the following scenario as a possible mech- 
anism of the anomalous temperature dependence of x(0) observed by Bobroff et. al.: (1) 



When the impurity scattering is in the unitarity hmit, the low-energy bound state appears 
around the impurity in d-wave superconductivity. (2) Then the AF spin fluctuations are 
enhanced around the impurity by the DOS originating from the bound state. (3) This local 
enhancement affects x(0) through the mode-mode coupling effect mentioned in the impurity 
effect-(B). As a result, x(0) increases with decreasing temperature around the impurity even 
in the superconducting state. 

In this paper, we investigate the possibility of this scenario in high-Tc cuprates. For this 
purpose, we consider a model two-dimensional (i^2_j^2-wave superconductor with strong AF 
spin fluctuations and with a nonmagnetic impurity. In this model, we calculate the static 
spin susceptibility around the impurity. We show that the local enhancement of x(0) really 
occurs around the impurity owing to the two effects, (A) and (B), noted in the above. 

This paper is organized as follows: In §2, we present our formulation. Nonmagnetic 
impurity effects on susceptibility are studied in §3, which is followed by summary in §4. 

II. FORMULATION 

A. Model two-dimensional d^2 _y2-waLve superconductor 

We consider an x A^-two-dimensional square lattice with a nonmagnetic impurity at 
Rimp = (0,0). We impose the periodic boundary condition in the x- and ?/-direction and 
put the lattice constant unity. The Hamiltonian is given by 

^ = I] [tij4aCja + h.c] - /i ^ n, + ^ Ui^Uii ^9^2 ^i^^jl + ^imp, (2.1) 

i,j,cr i i {ij) 

where Qo- is the annihilation operator of an electron with a-spin at the i-th site, Uia- = cj^Cia, 
and rii = J^a ^ia'-, U and g represent the on-site Coulomb repulsion and the nearest-neighbor 
pairing interaction, respectively. In the paring term in eq. ( |2.1|) , {ij) means the summation 
over the nearest-neighbor sites. 

In the hopping term in eq.(^TT]), we put tij = —t = —1 for the nearest-neighbor hopping 
and tij = t' = 0.25 for the next nearest-neighbor ones. (In this paper, t = 1 is the unit 
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of energy.) In addition, the chemical potential fi is chosen so that the number density of 
electrons can he n = 0.95 at sites far away from the impurity. Then we obtain a nearly 
antiferromagnetically nested Fermi surface, which gives strong AF spin fluctuations when 
f/>0. i- 

The last term in eq. ( |2.1| ) describes the nonmagnetic impurity scattering. A recent exper- 
imental analysis clarified that the suppression of Tc by Li-doping in YBa2Cu306+a; is almost 
equal to the case of Zn-doping. |]l5l Since Zn is considered a unitarity scatterer, this result 



indicates that Li should be also regarded as this type of scatterer. For this reason, we put 
Vi = 100 as the strength of the impurity potential at Rimp- In addition, we assume that the 
impurity potential at the nearest-neighbor sites of the impurity is finite in order to examine 



the effect of the extended impurity potential in real space. |T8| The resulting i^imp is then 
given by 

H,n.p = Vinn,^^ + V2 (2.2) 

r,=(0,±l),{±l,0) 

where is the position vector of the j-th site. 

We apply the mean field approximation to the pairing term in eq. Then this term 

is reduced to 

-gY. ^i^^Jl Jll^iAAi + (2-3) 

where A^- = —g{cjic^) is the superconducting order parameter. In the present case 
(t = = 0.25, n = 0.95), the most stable superconducting state is the dx^^-y^-wme su- 
perconductivity. The order parameter then satisfies A^j = Ajj and Ajj+j^ = — Aj.j+j, in the 
region far away from the impurity, {i + x{y) is the nearest-neighbor site of the i-th site in 
the x(?/)-direction.) 

B. Spin susceptibility 



We calculate the following static spin susceptibility within the Hartree-Fock random 
phase approximation (HF-RPA): 
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X(r„r,)=i/ dte-'\[a^it),a!m), (2.4) 
Jo 

where cr|(t) = — riiiit) is the spin operator in the Heisenberg representation. In 

calculating the susceptibility, we extend our previous method [§] to the case of the super- 
conducting state: (1) We replace the Coulomb interaction term in eq. ( |2.1| ) by the one under 
HF approximation: 

U E ^i^^n ^ J T.(^^)^i - I E(^')^r + T.[^4Ai + (2-5) 

i i i i 

where we have dropped constant terms. In eq. Ai = U (ci|Cj|) is the s-wave supercon- 

ducting order parameter which is induced around the impurity by the Coulomb interaction 
and the d^2_y2-wave superconductivity. ||T9|^T| (2) Then the Hamiltonian can be diagonal- 
ized by the Bogoliubov transformation in real space: (cjp clp ■ ■ ■ , cjy2|, ci|, C2|, ■ • ■ , 0^2^) = 
ilh 72) ■ ■ ■ 5 72Ar2)W^^ where ly is a 2N^ x 2A^^-unitary matrix while 7^ represents the anni- 
hilation operator of an eigen-state with an energy Ei. (3) The spin susceptibility under HF 
approximation (= x^) is calculated from 



(2.6) 

In eq. ([^.6D, f{E) is the Fermi distribution function. (4) The spin susceptibility under 
HF-RPA is given by x = [1 - (t//2)xT^X°> where x = {x(ri,rj)} and x° = {x°(ri,rj)} 
are N"^ x A^^-matrices. (5) We execute the Fourier transformation in terms of the relative 
coordinate, — r^, in x(rj, r_,): 

x(q,R) =Ee'''''^''"'"^^x(2R-r,,r,-), (2.7) 
j 

where R = (rj + rj)/2, and the summation is taken over all lattice sites. Since R is not 
always at a lattice site but it can be located between sites, we define the "on-site" spin 
susceptibility by M 



x(q,R) = x(q,R) + ^ E x(q,R + d) 

<i={0,±l/2), (±1/2,0) 

+ - E x(q,R + d). (2.8) 

d=(l/2,±l/2),(-l/2,±l/2) 

We use eg. (p. 81) in the following analyses. 

In numerical calculations, we determine Aj^, Aj, (nj), and {a^) self-consistently. In our 
previous paper, a large lattice (A^ < 52) was necessary in order to avoid finite size effects 
in calculating x the normal state. On the other hand, this effect was found weak below 
Tc, so that we put = 30. The A^-dependence of x is then very weak except for the case of 
g = 0.9 in Fig. 3: In this case, weak size dependence still remains just below Tc; however, 
it does not change our conclusions. As for x above Tc, we used the results in our previous 
paper, g 

III. RESULTS 

A. Local enhancement of spin susceptibility below Tc 

Figure 1 shows the temperature dependence of x(q; R-) in the case of V2 = 10. In 
Fig. 1(a), x(0, R = (1,1)) decreases just below Tc = 0.398 owing to the suppression by 
superconductivity and then increases at lower temperatures. The increase is more remarkable 
than that in the case when the superconductivity is absent. (Compare x(0, R = (1, 1)) with 
the result, " (l,l):normal" : The latter is obtained by putting g = 0. ||2^) The temperature 
dependence of x(0, R = (1, 1)) qualitatively agrees with the Knight shift experiment on Li 

in YBa2(Cui_,Li,)3 06.97. i 

On the other hand, the temperature dependence of x(0,R = (5,5)) is almost the same 
as that in the absence of the impurity (solid circles). This means that the anomalous 
temperature dependence of x(0, R) occurs in the vicinity of the impurity only. Indeed, 
when we examine the spatial dependence of x(0, R), we find that x(0, R) is enhanced only 
near the impurity site as shown in Fig. 2(a). 
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We also obtain an anomalous behavior in the AF susceptibility. Below Tc, Fig. 1(b) 
shows that x(Q) R- = (1; 1)) (Q = (^j ^)) is larger than that in the absence of the impurity 
(solid circle) and increases remarkably at low temperatures. This anomalous temperature 
dependence is not obtained at R = (5,5). At this site, x(Q)R') is almost temperature- 
independent below Tc. Namely, as in the case of the uniform susceptibility, the enhancement 
of the AF susceptibility occurs only in the vicinity of the impurity (Fig. 2(b)). 

Figure 3 shows the effect of the strength of the pairing interaction on the anomalous 
increase of x(0,R = (1,1)). In this figure, we find that the temperature region where 
x(0, R) is suppressed by superconductivity becomes narrower as g becomes smaller. In 
particular, at g = 0.9, x(0,R = (1,1)) continues to increase with decreasing temperature 
below Tc = 0.11 as if the system is not in the superconducting state. Since the AF spin 
fluctuations are stronger at lower temperatures. Fig. 3 indicates that the local enhancement 
of the uniform susceptibility is more remarkable when the AF spin fluctuations are stronger. 

Experimentally, |1| the Knight shift on Li decreases just below Tc and then increases 
at lower temperatures in YBa2(Cui_^Lij.)3 06.97 (Tc = 79.5[K]). On the other hand, the 
Knight shift continues to increase below Tc in YBa2(Cui_^Li^)306.6 (Tc = 41 [K]). Since the 
AF spin fluctuations are considered stronger in the latter, our result is consistent with this 
experimental tendency. 

B. Effects of low-energy bound state and mode-mode coupling 

In order to explain the mechanism of the anomalous temperature dependence of the spin 
susceptibility, we briefly review the mode-mode coupling effect which was discussed in our 
previous paper: When we consider Xij up to 0{U) in a continuum system, we obtain ^ 

X(q, R) = X°(q, R) + |e^[^''^-^«i-^''i-^«2l/(qi, Ri)x°(q2, R2), (3.1) 

where we put qi q and Rj — > R (i = 1,2) after the gradients are executed. When the 
system is homogeneous, V-rX° = 0- Then the gradient in terms of q does not work in eq. 
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In this case, is determined by x*^ at Q only. On the other hand, when VhX° 0, 
V,X°(ci) affects which means that with q' 7^ q is necessary in obtaining x(q) 

(mode-mode couphng). 

Within the neglect of the mode-mode coupling effect, the susceptibility is given by 

Figure 4 shows the comparison of R) with x(q, R) which is given by eq. (|2.8| ) with % 
being replaced by x'- In Fig- 4(b), we find that the enhancement of the AF susceptibility 
at 0.2 < T < Tc can be reproduced by x apart from small deviation. When we calculate the 
local DOS which is defined by 

iV(^,r,) = J2 mm\'S{u;-Ej, (3.3) 

m=l 

we obtain large intensity around co' = OatR = (1,1) originating from the low-energy bound 
state as shown in Fig. 5. [^,|| Since x'{^^ R-) is affected by the local DOS at R, the low-energy 
bound state is considered the origin of the local enhancement of the AF susceptibility. 

On the other hand, the local enhancement of the uniform susceptibility cannot be ex- 
plained within the neglect of the mode-mode coupling (Fig. 4(a)). Thus, as in the case of the 
normal state, |^ the mode-mode coupling is crucial in the case of the uniform susceptibility. 

Roughly speaking, when spin fluctuations are enhanced locally in real space, this phe- 
nomenon is described as the enhancement of the spin susceptibility at various q in momen- 
tum space. In the present case, although mainly the AF spin susceptibility is enhanced by 
the effect of the modification of the local DOS around the impurity, the uniform suscepti- 
bility is also enhanced because of the fact that the enhancement of the spin fluctuations is 
localized in real space. 

Figure 6 shows the temperature dependence of the spin susceptibility in the case of the 
5-functional impurity potential (V2 = 0). In this case, although both x(0, R) and x(Q,R) 
at R = (0, 0) are suppressed above Tc compared with those at R = (5, 5), they increase far 
below Tc as in the case of V2 = 10. At low temperatures, they are enhanced locally around 
the impurity (Fig. 7). 
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In the normal state, the spatial dependence of the impurity potential is crucial in obtain- 
ing the Curie- Weiss-like temperature dependence of x(0) around the impurity. 0,|1^,|I^ For 
example, the anomalous temperature dependence is not obtained when V2 = 0: In this case, 
the local enhancement of DOS around the Fermi level does not occur in contrast to the case 
of V2 7^ 0. Then the AF spin fluctuations are not enhanced around the impurity, so that 
the anomaly in x(0) is also absent. On the other hand, in the dx^^y^-wave superconducting 
state, since the low-energy bound state exists around a unitarity scatterer even in the case 
of V2 = 0, the local DOS has large intensity around u = near the impurity (inset in Fig.6). 
Then the AF spin fluctuations become stronger using this large intensity of the local DOS, 
so that the uniform susceptibility is also enhanced by the mode- mode coupling effect. 

Comparing Figs. 1 and 2 (V2 = 10) with Figs. 6 and 7 (V2 = 0), we find that the 
local enhancement of the spin susceptibility is more remarkable in the former case. Namely, 
detailed spatial dependence of the impurity potential is crucial in considering the local 
enhancement of the spin susceptibility quantitatively. 



C. Local moment 



In Fig. 1, the AF phase transition occurs at T = Tm — 0.1. Below Tm, x(Q)R') 
R = (1, 1) decreases with decreasing temperature, and the spontaneous magnetic moment 
{{(Jz)) grows at this site as shown in Fig. 8. This antiferromagnetism is inhomogeneous 
and is localized around the impurity (Fig. 9); the induced moment is largest at (±1,±1). 
Because of this inhomogeneity, we obtain a finite total magnetic moment when we sum up 
(cr^) spatially. 

Since the present calculation is based on the mean field approximation, this local mag- 
netic instability itself is an artifact. We should regard Tm as a characteristic temperature 
at which the local moment is formed around the impurity. However, when the impurity 
concentration is finite, it might be possible that the local antiferromagnetism around each 
impurity connects to each other three-dimensionally and then the magnetism is stabilized. 
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Physical properties below Tm and the possibility of the stabilization of this antiferromag- 
netism remain as our future problems. 

IV. SUMMARY 

In this paper, we investigated the uniform susceptibility and the AF one around the non- 
magnetic impurity in (i^2_j,2-wave superconductivity. We calculated the spin susceptibihty 
in real space taking into account the effects of the extended impurity potential, dx2_y2-w8ive 
superconductivity and the Coulomb interaction. 

When the impurity scattering is in the unitarity limit, the low-energy bound state ap- 
pears around the impurity, which increases the intensity of the local DOS around u = 0. 
This increase enhances the AF spin fluctuations around the impurity at low temperatures. 
This phenomenon affects x(0) thought the mode-mode coupling effect, so that x(0) also 
increases at low temperatures. In addition, the enhancement of x(0) is more remarkable 
when the AF spin fluctuations are stronger. In particular, when the AF spin fluctuations 
are strong to some extent, x(0) near the impurity continues to increase below Tc without 
being suppressed by superconductivity. These results can explain the recent Knight shift 
experiment in YBa2(Cui_a;Li^)306+y. 

With regard to the effect of the extended impurity potential in real space, we showed 
that, apart from quantitative difference, the local enhancement of x(0) is also obtained when 
the impurity potential is finite only at the impurity site. This result is in contrast to that 
in the normal state; above Tc, the local enhancement of x(0) is not obtained as far as we 
use the ^-functional impurity potential. The difference originates from the existence of the 
low-energy bound state in the case of the (ia;2_y2-superconducting state; it enhances the AF 
spin fluctuations, which increases x(0) locally. 

We also showed within the mean field theory that the local antiferromagnetism may 
appear around the impurity below Tm. Although this antiferromagnetism itself is an artifact 
of the present mean field approximation, we can still expect a formation of a local moment 
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around Tm- Since the present mean field theory is inadequate to examined the region below 
Tm, a more sophisticated theory is necessary in order to clarify magnetic properties below 
Tm- This problem remains as our future problem. 
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Figure Captions 

Fig.l: Temperature dependence of spin susceptibility in the case of V2 — 10. (a) uniform 
susceptibility, x(0,R). (b) AF susceptibility, x(Q,R). Here, for example, (1, 1) is the 
result for R = (1, 1). In this figure and in Figs. 3, 4, and 6, 'normal' means that the 
system is in the normal state {g — 0), while 'pure' is the case in the absence of the 
impurity. We put U — 2, g — 2 and Vi — 100; in this case, we obtain Tc — 0.398. 
These values arc used also in the following figures except for Fig. 3. At T = Tm ~ 0.1, 
the AF phase transition occurs. {x{<l, R) below Tm is shown in Fig. 8.) 

Fig.2: Spatial dependence of x(q, R) at T = 0.1116 in the case of V2 = 10. 

Fig.3: (^-dependence of the uniform susceptibility at R = (1, 1). We put V2 = 10 and U = 2. 
Transition temperatures are = 0.398 = 2), Tc = 0.21 {g = 1.3), and Tc = 0.11 
U = 0.9). 

Fig.4: Comparison of x(q, R) with x(q, R) at R = (1, 1). The cusp in x(q, R) at T ~ 0.1 is 
due to the AF phase transition at Tm. 

Fig.5: Local DOS at T = 0.15. In calculating the DOS, we put N — 50 and added small 
imaginary part T = 0.2 to eigen-energies. 

Fig. 6: Temperature dependence of spin susceptibility around the impurity in the case of 
V2 = 0. The inset shows the local DOS at T = 0.15. 

Fig. 7: Spatial dependence of x(q, R) at T = 0.05 in the case of V2 = 0. 

Fig. 8: Temperature dependence of x(Q,R) above and below Tm at R = (1,1) {V2 = 10). 
The inset shows the spontaneous magnetic moment (cr^) at R = (1, 1). 

Fig. 9: Spatial dependence of (a) the d^^-y^-waye order parameter and (b) the spontaneous 
magnetic moment. We put V2 = 10 and T — 0.05. In this figure, Ad(rj) = 
i E«i=(±i,o),(o,±i) K„r,+a ■ i-l)'^ and a^{r,) = af ■ e-I^^+^^l (r, = (x,, y,))- 
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